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I. INTRODUCTION 

The analysis of the J/ip dissociation cross section is important for understanding the suppression of J/ip 
production observed in Pb-Pb collisions by the NA50 Collaboration at the CERN Super Proton Synchrotron 
(SPS) [1]. There are a number of theoretical calculations on the cc + light hadron cross sections (see, e.g., the 
review Ref. [2]). However, they give widely divergent results, which implies that one is still far away from a 
real understanding of the scattering mechanism. The nonrelativistic quark model has been applied in [3] and 
[4-6] for the calculation of the cross sections for the dissociation processes cc + qq — > cq + qc. The calculated 
cross sections for the reactions J/ip + n — > DD, D*D, DD*, D*D* have the following common features: they 
rise very fast from zero at threshold to a maximum value and finally fall off due the Gaussian form of the 
potential. The magnitude of the maximum total cross section was found to be w 7 mb at \fs w 4.1 GeV in [3] 
and a somewhat smaller value of w 1.4 mb at yfs w 3.9 GeV in [4-6]. 

Another approach to studying the charm dissociation process started with the model proposed by Matinian 
and Miiller [7] . They assumed that the dissociation cross section is dominated, in the t channel, by the D meson 
exchange. A generalization of this approach can be found in [8-13] where an effective chiral SU(4) Lagrangian 
was employed. Such an approach seems to us quite problematic for the following reasons: (i) SU(4) is a badly 
broken symmetry, and (ii) some of the couplings in the chiral SU(4) Lagrangian are unknown. Nevertheless, 
this is a relativistic approach that allows one to study the above processes in a systematic fashion. In this 
framework, the dissociation cross section of J/ip by light hadrons is predicted to be, near the physical threshold, 
in the range 1-10 mb. Moreover, it is interesting to refer to the paper [12], where the J/ip dissociation by 7r 
and p mesons was examined in the meson exchange model [12] and compared to the quark interchange model. 
The authors of this paper found that the meson exchange model could give predictions similar to those of the 
quark interchange models, not only for the magnitudes but also for the energy dependence of the low-energy 
dissociation cross sections of the J/ip by it and p mesons. 

It appears that the microscopic quark nature of hadrons is important in the charm dissociation processes. 
The first step is to calculate the relevant form factors corresponding to the triple and quartic meson vertices in 
the kincmatical region of the dissociation reaction. QCD sum rules have been used in Refs. [14-17] to evaluate 
those form factors and to determine the charm cross section. The cross section was found to be about 1 mb at 
y/s w 4.1 GeV with a monotonic growth when the energy is increased [17]. 

An approach based on the 1/N expansion in QCD combined with the Regge theory gives, for the total 
dissociation cross section, a value of a few millibarns near to y/s w 4 GeV [18]. 

We also mention the work of Deandrea et al, where the strong couplings J/ipD*D* and J/ipD*D*w were 
evaluated in the constituent quark model [19]. Finally, an extension of the finite-temperature Dyson-Schwinger 
equation approach to heavy mesons and its application to the reaction J/ip + it — > D + D was considered in 
[20]. 

We employ a relativistic quark model [21] to calculate the charm dissociation amplitudes and cross sections. 
This model is based on an effective Lagrangian which describes the coupling of hadrons H to their constituent 
quarks. The coupling strength is determined by the compositeness condition Zh = [22] where Zh is the 
wave function renormalization constant of the hadron H. One starts with an effective Lagrangian written 
down in terms of quark and hadron fields. Then, by using Feynman rules, the S-matrix elements describing the 
hadronic interactions are given in terms of a set of quark diagrams. In particular, the compositeness condition 
enables one to avoid a double counting of the hadronic degrees of freedom. The approach is self-consistent 
and universally applicable. All calculations of physical observables are straightforward. The model has only a 
small set of adjustable parameters given by the values of the constituent quark masses and the scale parameters 
that define the size of the distribution of the constituent quarks inside a given hadron. The values of all fit 
parameters are within the window of expectations. 

The shape of the vertex functions and the quark propagators can in principle be found from an analysis of 
the Bethe-Salpeter and Dyson-Schwinger equations as was done, e.g., in [23]. In this paper, however, we choose 
a phcnomcnological approach where the vertex functions are modeled by a Gaussian form, the size parameter 
of which is determined by a fit to the leptonic and radiative decays of the lowest- lying light, charm, and bottom 
mesons. For the quark propagators we use a local representation. 

We calculate the amplitudes and the cross sections of the charm dissociation processes 

J/ip + 7T -> D + D, 
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J/V' + tt -» £>* + £> (D*+D), 
J/V' + tt -» D* +W. 

These processes are described by both box and resonance diagrams which can be calculated straightforwardly 
in our approach. We compare our results with the results of other studies. 

The layout of the paper is as follows. In Sec. II we briefly discuss our relativistic quark model. In Sec. Ill 
we outline the calculational technique of the arbitrary n-point one-loop diagrams with local propagators and 
Gaussian vertex functions. We give explicit result for the triangle and box diagrams which are the building 
blocks of the charm dissociation amplitudes. In Sec. IV we calculate the charm dissociation amplitudes and the 
total cross sections. In Sec V we perform the numerical analysis and give our prediction for the cross sections. 



II. THE MODEL 



The coupling of a meson H to its constituent quarks qi and q 2 is determined by the Lagrangian 

£fnt ( x ) = 9hH{x) JdxiJdx 2 F H (x,x\,x 2 )q 2 {x2)T H \ H q 1 {xi) + h.c. . (1) 

Here, Xh and T H are Gell-Mann and Dirac matrices which describe the flavor and spin quantum numbers 
of the meson field H(x). The function Fjj is related to the scalar part of the Bethe-Salpeter amplitude and 
characterizes the finite size of the meson. To satisfy translational invariance, the function Fh has to satisfy the 
identity Fh(x + a, x\ + a, x 2 + a) = Fh{x, x\,x 2 ) for any four-vector a. In the following we use a particular 
form for the vertex function, 

F H (x,x 1 ,x 2 ) = S(x - c\ 2 x 1 - c\ 2 x 2 )§ H ((x 1 - x 2 f) , (2) 

where <&h is the correlation function of two constituent quarks with masses mi, m 2 and c\- — mk/{mi + rtij). 

The coupling constant gn in Eq. (1) is determined by the so-called compositeness condition originally pro- 
posed in [22] and extensively used in [24]. The compositeness condition requires that the renormalization 
constant of the elementary meson field H(x) is set to zero, 

Z * = !" %&h{M%) = 0, (3) 
where Y1' H is the derivative of the meson mass operator. In order to clarify the physical meaning of this 

1 /2 

condition, we note that Zjj is also interpreted as the matrix element between a physical particle state and the 
corresponding bare state. For Zh = it then follows that the physical state does not contain the bare one and 
is described as a bound state. The interaction Lagrangian in Eq. (1) and the corresponding free Lagrangian 
describe both the constituents (quarks) and the physical particles (hadrons) which are bound states of the 
quarks. As a result of the interaction, the physical particle is dressed, i.e., its mass and wave function have 
to be renormalizcd. The condition Zh = also effectively excludes the constituent degrees of freedom from 
the physical space and thereby guarantees that a double counting of physical observables is avoided. The 
constituent quarks exist in virtual states only. One of the corollaries of the compositeness condition is the 
absence of a direct interaction of the dressed charged particle with the electromagnetic field. Taking into 
account both the tree-level diagram and the diagrams with the self-energy insertions into the external legs 
yields a common factor Zh which is equal to zero. We refer the interested reader to our previous papers 
[21, 24, 25] where these points are discussed in more detail. 

We briefly discuss the introduction of the electromagnetic field into the nonlocal Lagrangian in Eq. (1) in a 
gauge invariant manner. This can be accomplished by using the path exponential [26] 

£?nt +el » = 9hH(x) Jdx^dx^Hix^^xM^e^^^^rHXHe-^'^^Ui^i), (4) 



where 
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P) = Jdz^(z 



(5) 



and where z M is a coordinate point on the path P. At first sight it appears that the results depend on the path 
P which connects the end points x and x\ in the path integral in Eq. (5). However, we need to know only 
derivatives of such integrals for the perturbative calculations here. Therefore, we use a formalism [26] which is 
based on the path-independent definition of the derivative of I(x, y, P): 



d 

Km dx" — I{x,y,P) = lim [I(x + dx, y, P') - I(x, y, P)\ 

dxt^O OX^ cfa;f-+0 



(6) 



where the path P 1 is obtained from P by shifting the end point x by dx. The use of the definition in Eq. (6) 
leads to the key rule 



d 
~dxM 



I{x,y,P) =A ll {x) , 



(7) 



which in turn states that the derivative of the path integral I(x, y, P) does not depend on the path P originally 
used in the definition. This allows one to construct the perturbation theory in a consistent way and guarantees 
the implementation of charge conservation and the Ward identities. 

As an example, we consider the transition J/ip — > 7 which is now described by two diagrams in Fig. 1: (a) 
the standard "bubble" and (b) the "tadpole" ones. The total matrix element has a manifestly gauge invariant 
form 



M^(P) = OTP 2 



(8) 



where 



3g.j/->p 1 



y a = 1 + 



4tt 2 
t 



^2 j dt (rTtp / ' da {y a ®j/A z a)-yb®j/*p( z b)} , 





, p 2 l + 3t/2-2t 3 a(l -a) 



z a = t[m 2 c -a(l-a)p 2 ] 



l + t\2 



-a) p 



Vb 

Zb 



2 2 
a p 



l + t 



(l + tf 4 
ap 2 



For the pseudoscalar and vector mesons treated in this paper the derivatives of the mass operators are written 



as 



%{P 2 ) = 
K(P 2 ) = 



J_ a J_ fd 4 k - 
2p 2P dp a J 4w 2 i 

1 



$|,(-r)tr 



1 5 S 1 (fi + c\ 2 ^ 5 S 2 ^-c 2 12 p 1 ) 



(9) 



p^p u 



j_ a j_ rd A k - 

2p 2? dp a J 4n 2 i 



:t> 2 v (-k 2 )tr 



7^i(^ + 4 ^S 2 m-c 2 2 f) 





(a) 



(b) 



FIG. 1: The diagrams corresponding to the J/ip — > 7 transition: (a) the usual one and (b) the "tadpole". 



The leptonic decay constants fp and fv are 
3g p f d A k ~ 



3flv 
4tt 2 



4tt 2 
d A k ~ , 



= fpp" 



k 2 )tr 



C12 f) 



my fv f-v ■ 



We use free fermion propagators for the valence quarks 
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(10) 

rrii- ft 

with an effective constituent quark mass rm. As discussed in [21, 24] we assume for the meson mass Mh that 

M H <m 1 +m 2 (11) 

in order to avoid the appearance of imaginary parts in the physical amplitudes. This holds true for the light 
pseudoscalar mesons but is no longer true for the light vector mesons. We shall therefore employ identical 
masses for the pseudoscalar mesons and the vector mesons in our matrix element calculations but use physical 
masses in the phase space calculation. This is quite a reliable approximation for the heavy vector mesons, e.g. 
D* and B*, where the hyperfine splitting between the D* and D and the B* and B, respectively, is quite small. 

The shape of the vertex functions and the quark propagators can in principle be found from an analysis of the 
Bcthc-Salpeter and Dyson-Schwingcr equations as was done, e.g., in [23, 27]. In this paper, however, we choose 
a phenomenological approach where the vertex functions are modeled by a Gaussian form, the size parameter 
of which is determined by a fit to the leptonic and radiative decays of the lowest- lying light, charm, and bottom 
mesons (see Table I) . Our previous studies of phenomena involving the low- lying hadrons have shown that this 
approximation is successful and reliable [21, 24]. We employ a Gaussian for the vertex function of the form 
<f>ij(k%) = exp(— k 2 E /K 2 H ), where fee is a Euclidean momentum. The size parameters K 2 H are determined by a 
fit to experimental data, when available, or to lattice results for the leptonic decay constants fp and fv where 
P = 7r, D, B and V = J/ip, T. Here we improve the fit by using the MINUIT code in a least-squares fit. The 
values of the fit parameters are displayed in Table II. The quality of the fit may be assessed from the entries 
in Table I. 



III. STRONG TRIANGLE AND BOX DIAGRAMS 



Transition matrix elements involving composite hadrons are specified in the model by the appropriate quark 
diagram. Here, we give explicit expressions for the integrals corresponding to the strong triangle and box 
diagrams shown in Fig. 2. 
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This model 


Expt. or lattice 




This model 


Expt. or lattice 


u 


130.7 


130.7 ±0.1 ±0.36 




0.272 GeV -1 


0.273 GeV -1 




159.8 


159.8 ± 1.4 ±0.44 


f.T/4> 


405 


405 ± 17 


Id 


211 


203 ± 14 
226 ± 15 


fa 


182 


173 ± 23 
198 ± 30 


fo e 


244 


230 ± 14 
250 ± 30 


f Bs 


209 


200 ± 20 
230 ± 30 




360 


360 


fr 


710 


710 ± 37 



TABLE I: The physical quantities used in least-squares fitting our parameters. The values are taken from the Particle 
Data Group [28] or from lattice simulations [29]. The value of /s c is our average of QCD sum rules calculations [30]. 
All numbers are given in MeV except for g 1T a 11 . 



Quark masses 


(GeV) 


Ah 


(GeV) 


Ah 


(GeV) 


m u = m d 


0.223 


K 


1.074 


As c 


1.959 


m s 


0.356 


Ax 


1.514 


Aj/y, 


2.622 


m c 


1.707 


A D = A Ds 


1.844 


At 


3.965 


m b 


5.121 


A s = A Ss 


1.887 







TABLE II: The fitted values of the model parameters. 



First, we will make the transformations which are common for all one- loop diagrams with local propagators 
and Gaussian vertex functions. The Feynman integral corresponding to a one-loop diagram with n propagators 
and, respectively, n vertex functions may be written in Minkowskii space as 



I„(Pl, ...,p n ) 



d 4 k 
47T 2 ' 



■tr 



Y[$i{-(k + v i+n ) 2 ) FiSii^ 



(12) 



where the vectors Vi are linear combinations of the external momenta Pi to be specified later on, k is the loop 
momentum, and Ti are Dirac matrices for the i meson (cf. Eq. (1)). The external momenta are all chosen as 

n 

ingoing such that one has Pi = 0. 

i=l 

The propagators can be written as 



$(/f+ ft) = (rm+ fi+ fa) ■ /dAe- ft(m? - (fc+, " )a) . 



(13) 



For the vertex functions one writes 

^(-(k + V t + n f) = e 0> + n-(k+V i+ nf 



1, ...,n, 



(14) 



where the parameters (3i +n = s, = 1/A? are the size parameters. One can then easily perform the integration 
over k: 



/<i 4 fc f 2 ™ 1 1 I 1 

U=l J f I P l<i<j<2n 



(15) 
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fc + c 



13 Pi J 



[k + pi + c^p 3 ] 



P4 



k +pi 



e 



fc — P2 — P3 



»»3 



fc — C 



12P2J 



* - P2 - cl 3 p 3 ] 



P3 



fe - P2 



0^. 



[fc + cj 4 pi] 



[fe - c; 2 p 2 i 



FIG. 2: The strong triangle and box diagrams. 



Here, j3 = ^2 Pi- The numerator vrn + }i + fa can be replaced by a differential operator in the following manner: 



fd 4 k -A- 

/ — oT TT Ti (mi+ }£+ fa)e 



0k 2 +2kr 



JJr 4 ( m,+ $ + £ # 

1=1 



i \ J_ e _ rv/3 



2 r ' y /? 2 



(16) 
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1 " / 11 

= ^ e " r2//3 n r ^ m »+^-^^+2 & 



We thus have 



o 



l<i<3<2n 



fir, ^+^-1/+!^ 



(17) 



»=i 
1 

x -tr 

4 



In 



where r = ^ /3j Vi . Finally, we effect some further transformations on the integration variables to get the 

i=i 

integral into a form suitable for numerical evaluation: 

• we use the formula 

n °° °° n °° / n \ 

JJ d/3 i f(f3 1 ,...,/3 n )= dte-'H da t S \f(t ai ,...,ta n ); 
{ i=i{ \ i=i / 

n n 

• we scale the t variable by t — > w i with to = /3 i+ „ = ^ s« and introduce the new variable /3j = 

i=l i=l 

(i = n + 1, 2n). 



We have 



* (T+7)2 II y d«i j ( i - ^ 

1=1 ^ 1=1 J 

' n / I I 

nr^+A- — f+- 



x — tr 

4 



where 



(18) 



f = t ^ a i V i + & Vl ■ 



i=l 



The z form in the exponential function is written as 



(19) 



^loc 



^2 = 



i—1 l<i<j<n 
n 2 n 

z— 1 j—n+1 1 <n 

X! & ft a v ■ 

n+l<i<j<2 n 



The matrix a,ij = (w, — vj) 2 depends on the invariant kincmatical variables. Explicit expressions for this matrix 
for the triangle and box diagrams are given in Tables III and IV, respectively. We will introduce the variable 
v = t/(l + t) (0 < v < 1) in Eq. (18) in what follows. 
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ffli2 = pi 


Ol3 = Pi 


123 = Pi 


an = (cj 3 ) 2 p? 


Ol5 = (cj 2 ) 2 pi 


oie = ci 3 p 2 - ci 3 c 23 pi + c 23 pi 


an = C? 3 pi - c} 3 C? 3 p? + c\ 3 pi 


/ 2 \ 2 2 
125 = {C 12 ) V'l 


/ 2 \ 2 2 
126 = (C 23 ) p :i 


134 = (c? 3 ) 2 pf 


135 = c? 2 pj - c\ 2 ci 2 pi + cj 2 pi 


136 = (c 23 ) 2 pi 


145 = C13 (Ci 3 - c\ 2 )pl 

+c\ 2 (ci 2 -cj 3 )pi + ci 2 c ; 3 pi 


146 = C13 (C13 — C 23 ) Pi 

+cl 3 (cis - c? 3 ) pi + c? 3 ci 3 pi 


2/2 2 \ 2 
156 = C12 (C12 - C 23 )P2 

+ci 3 (c| 3 - c? 2 ) pi + c? 2 ci 3 p? 



TABLE III: The matrix a,ij = (uj — Vj) 2 for the triangle diagram. Here, vi = 0, v 2 = — P2, «3 = pi, «4 = c 13 pi, 

f 5 = ~c\ 2 Pi, V 6 = Pl + C 23 P3- 



2 

112 = P2 


ai3 = (p2 +P3) 2 


ai4 = Pi 


123 = Pi 


124 = (pi +P2) 2 


134 = P4 


1 1 \ 2 2 
115 = (,Ci 4 J Pl 


1 1 \ 2 2 

ai6 = (c 12 J p 2 


an = C23 pi — ci 3 C23 pi 

+c| 3 (P2 +P3) 2 


us = C34 pi - C34 C34 pl 

+44 (Pl +P4) 2 


4 2 14 2 
125 = C14P2 ~ C14 C14P1 

+c} 4 (pi +P2) 2 


1 2 \ 2 2 
126 = (C12J p 2 


( 2 \2 2 
127 = (C 23 J p 3 


4 2 3 4 2 
128 = C34 p 3 — C 34 C34 p 4 

+C34 (P3 +P4) 2 


12 14 2 
135 = C14 p 4 — C14 C14 Pi 

+cf 4 (pi +P4) 2 


136 = cl 2 Pi - cj 2 C 2 2 pi + C? 2 (P2 + P'lf 


/ q \ 2 
137 = (C 23 j p 3 


138 = (C34J P4 


1 4 \2 2 
(145 = (C14J pi 


2 2 12 2 
a 4 6 = C12 Pl - C12 C12 p 2 

+c} 2 (pi +P2) 2 


147 = ci 3 pi - c| 3 ci 3 pi 

+Cm(P3+P4) 2 


/ 4 \2 2 
148 = (C34J P4 







TABLE IV: The matrix dij = (vi — Vj) 2 for the box diagram. Here, vi = 0, v 2 = —p 2 , V3 — —p 2 — P3, V4 — pi, 

W 5 = C 14 pi, V(j = -c} 2 P2, V 7 = ~Pl - ci 3 P3, V S = pi + c| 4 p 4 . 



Some further remarks are appropriate. One can see that the existence of the loop integral Eq. (18) is defined 
by the local a form z\ oc . Even if we apply the constraints Eq. (11) it does not guarantee that this form is 
always positive. For instance, in the simplest triangle diagram with p\ = p\ = p\ = p 2 and mi = m 2 = m 3 = m 
the form z\ oc = if p 2 = 3m 2 . Such singularities in the Feynman diagrams with local propagators are called 
anomalous thresholds (see the discussion in Ref. [31] and other references therein). The situation is much more 
complicated in the case of the box diagrams. Obviously, we cannot consider the annihilation processes when 
the energy is large enough to go beyond the normal thresholds corresponding to quark production. However, 
in the case of the dissociation processes considered here, the local form z\ oc is always positive, and we are able 
to make self-consistent and reliable predictions for physical observables. 
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D+ 




J/* 



D° J/* 




D° TV 



D+ J/9 





J/* 



D° J/ 




D° TV 



D*+ J I '* 





J/* 




D*° TV 



D*+ J/* 




FIG. 3: The Feynman diagrams describing the charm dissociation processes J/tp+-K + — » D + +D°, D* + +D°, D* + +D* °. 
IV. THE J/V» DISSOCIATION AMPLITUDES AND CROSS SECTIONS 



In our approach the dissociation processes J/ip + n + — > Z? + 13°, 13* + + D°, and D* + + D*° are described 
by the diagrams in Fig. 3. 

Our momentum labeling is defined by 



JMpi) + ir+(p 2 ) -> £>+(«!) + £>»(») 



(20) 



where D+ = D+ or D* + , D\ = D° or D*°, p\ = m 2 H1 = m 2 ^, V \ = m% 2 = ml, q\ = m 2 H3 = m 2 D+ {m 2 D , + ), 
q 2 = m 2 HA = rr2 m {m 2 WT ). 

The Mandelstam variables are defined in the standard form 

s = (pi +P2) 2 = (qi + 92) 2 , 
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t = (pi- qi) 2 = (P2 - qi) 2 -, 
u = (pi - qi) 2 = (P2 - qi) 2 , 



where s + t + u = m 2 H1 + m 2 H2 + m 2 H3 + m 2 Hi . 

Cross sections are calculated by using the formula 



°^=vLwrJ itm ' M2 (21) 

t_ 

where M(s, t) is an invariant amplitude and 

t± — (-El, cm — Es iCIn ) — (pi lC m T 9l,cm) 2 , 

_ a + m 2 Hl - m 2 H2 _ a + m 2 H3 - m 2 H4 

D l,cm n i— j -t^cm — n r~ ' 

2a/S 2a/S 

\ 1 / 2 (s,m 2 H1 ,m 2 H2 ) A 1//2 (s, m 2 H3 ,m'j I4 ) 
Pl.cm — t— = > 7i ~ — 



V- i yr,cm <-< r- 

The reaction threshold is equal to s = (mjf 3 + m^) 2 . Note that Eq. (21) contains the statistical factor 1/3 
which comes from averaging over the initial state J/ip polarizations. 

The dissociation processes are described by both the box and the resonance diagrams as shown in Fig. 3. 
The resonance diagrams depend explicitly only on the t or u variables whereas the box diagrams are functions 
of s and t. 

In the following three subsections we write explicitly the amplitudes for the processes J/ip + n + — > D + + D°, 
D* + + D°, D* + + D*° in terms of form factors; all the analytical expressions for them are reported. 

A. The channel J /if) + n+ — > D + + D° 

The invariant matrix element is written as 1 

M(VPPP) = M n (VPPP) + M res (VPPP), (22) 

where 

M n (VPPP) = e Ml (pi)-e" lPlPa9a Fvppp(s,t), (23) 
M res (VPPP) = e m (pi) ■ { - Ml (PVP) £># (P2 - 92) M^(WP) 
- M« b (PPV) (pa - gi ) M^f (VPV)} , 

M£ Q (PVP) = ^ ^pvp (*) > M^(WP) = e ^w(w-«) fWp(t), 

M£ b (PPV) = p« F^(u), M£f (VPV) = e fi/»Pi(w-90 F$ v (u). 

The common minus sign is due to the extra quark loop in the resonance diagrams as compared to the box 
diagram. This extra sign is very important, as will be seen later on, since the extra sign leads to a construc- 
tive interference of resonance and box graph contributions. The loop momenta directions come from explicit 



Note that the symbol e a ^ vq = e a f 3ll ' J p tl q v , where e a fiw ; s the Levi-Civita tensor. Moreover, the order of indices V and P of the 
form factors should be read looking at the Fig. 3 and following the arrows backward in the loops. 
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calculation of the S-matrix elements and are shown in Fig. 3. The vector (D*) and pseudoscalar (D) meson 
propagators are given by 



-g a ^ + p a p f3 /m 2 D , 



m 2 D , — p 2 



i 2 D D (p) 



-1 



m 2 D — p 2 



For the form factors one obtains 

fVppp(M) = 9j/,i>9d9* [d&n • cxp(~w ■ z n (s,t)) ■ w 2 (ra c (1 + d 2 ) - m q d 2 ), 



PVP 



(*) 



VVP 



(*) - 



PPV 



(«) 



4pv(«) 



5tt ffn* ffn / rfcr A • cxp(-wz A (t)) ■ [w {-m q m c - &i 
+m£ &2 ( &1 + ft 2 _ i) _ ^ b j &2 
+t& 2 (6i & 2 + &?-l + fo 2 )) 
-(l-v)(l + 3b 1 )}, 

9j/iP 9d* 9d Jda A • exp(-wz A (t)) ■ [w (m c (b 2 - 1) - m, & 2 )], 

9tt gD' gD JdaA ■ exp(-wzA(u)) ■ [to (m q m c + m 2 (1 - b\ - b 2 ) 

+ml (-2 &! & 2 + &! 6^ + &i + 2 6? b 2 - 2 b\ + &?) 
+m|> (-2 6i 6 2 + 2 6x &1 + 6? fe 2 + b 2 - 2 + b\) 
+u(2b 1 b 2 -b 1 b\-b\ b 2 )) 
+4 + 3 h v - 3 h + 3 6 2 v - 3 6 2 - 4 w] , 

9j/ip9DgD* / rfcr A • cxp(-toz A (u)) • [w (-m c + o 2 (m c - m,))]. 



The integration measures are defined by 



fa J '"(i-t 



JdaA = ^ Jdv (j^j J d*a6 (l - £>V 

V : I . 



The coefficients bi (triangle diagrams) and di (box diagrams) are given by 

h = va 3 + {l-v)(s- L c\ 3 + s 3 cl 3 ), 
b 2 = va 2 + (l-v)(s 2 c\ 2 + s 3 cl 3 ), 



(24) 



di = va 4 + (1 - v) (§i c{ 4 + S4C34), 

d 2 = -v(a 2 + a 3 ) - (1 - v) (§2 c\ 2 + s 3 + s 4 034), 

d 3 = -v a 3 - (1 - v) (s 3 c 23 + s 3 + s 4 Cg 4 ) . 

The values of za , 2n j and w are defined by the relevant diagrams (size parameters and quark and hadron 
masses). 
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B. The channel J /if) + tt + — > D* + + D° 

The invariant matrix element can be written as: 

M(VPPV) = M n (VPPV) + M res (VPPV) , (25) 

with 

M u (VPPV) = e w (pi)e Ma (gi) • (f>f pf 4ppv(^ *) + <Z2 1 pf 4ppv(s, *) 

+ tf 1 V? F^% v (s,t) + q^p^F^ PV (s,t)+g^F^ PV (s,t)) , 

M rcs (VPPV) = e^( Pl )eM 

x (- M£ (PVP) D a Jt( P2 q 2 ) M^(WV) - M^(PW) D#(P2 - gi) (VPV) 

+ M^ c (PVP)£> D (p2-?i)M^ c (VPP)) . 



The expressions for the form factors are given in the Appendix. 



C. The channel J/ip + 77+ — > D* + + D* 



The amplitude for the process J/ip + ir + — > D* + + D* can be written as 

M(VVPV) = M n (VVPV) + M ros (WPV), (26) 

where 

M n (VVPV) = e Ml (pi)e M2 (<Z2)e M 3(<Zi) 

x (jWm^bmi F« pv (s, t) + g^e?"*™ F$ PY (s, t) + gi*^!*™ F^ PY (s, t) 

+P^e PlPJ ' ,1 ' ,I 4vpv(M) +Pf e^" 1 " 3 4vpv(M) + pf 4vpv(M) 
+ ^ £ p lP ^^s F$ pv ( a ,t) +^£^» 4vpv(M) + g^ 1 ^ 2 " 3 4vpvM) 

fc -'VVPVv*' V ' "l fc r VVPVV s i L ) ' Pi fc r WPVV 4 i V 

+p ^ £P , m3 F^ PY (s, t) + j^^*"^ 4vpv(^ *) + g^ 292 ^ 3 4vpv( s - *) 

+e PiMiM a « 4^(5, *) + 4vP V (^ *)) • 



When taking into account parity invariance one knows from helicity counting that there are only 14 inde- 
pendent amplitudes in the (VVPV) case. The above set of 17 amplitudes are in fact not independent. They 
can be reduced to a set of 14 independent amplitudes by making use of the Schouten identity (see, e.g., [32]) 

+ cycl.(M 1 ,M2, Ms, M4, Ms) = 0. 

We have made use of the Schouten identity as an additional check on our numerical calculations. 

Moreover, the resonance term, M res (VVPV), and all the structures involved are reported in the following 
expressions: 

M res (VVPV) = e w (ft)e, 2 fe)^(9i) 

x (- MZ 2 (PW) D a £ (pa - q 2 ) M^ 3 (VVV) - MZ 3 (PW) (p 2 - gi) M£f ^ 2 (VVV) 
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+ M£ C (PPV) D D ( P2 - q 2 ) M£f (VPV) + M£ d (PVP) D D ( P2 - ?i) M£f (VVP)) , 
M^ Q 2 (PVV) = e Mp2 - q2)f * 2a F$ v {t) , 



^ W (WV) - (pa - 9 2 )^^pf i##(t) + (P2 - (P2 - ?2)M 3 4vv(*) 

+ .9 MlM P? ifvv (*) + (P2 - q 2 f F^ v (t) 
+ 4vv (*) + 5 M3/3 (P2 - ga)" 1 4vv W - 



M^ b 3 (PVV) = e^^-^^J^vCu), 

M£f"(VW) - (pa - gO^f 4vv(^) + (P2 - 9i) Ml (P2 - gi)V 2 4vv(«) 
+ .9 M1M2 P? ifvv («) + 5 mM2 (P2 ~ 9i)^ 



M£ C (PPV) = p^F^t), 
M^f(VPV) = £ piqi ^ 3 F$ y (t) , 

M£ d (PVP) = pf ^p(") 
M^f(WP) = e Pl<?2 ^^F v ^ P (u). 

All the expressions for the form factors are reported in the Appendix. 

Note that the J ftp dissociation amplitudes are not equal to zero when contracted with the four- momentum of 
the J/ip except in Eq. (22) which involves the Levi-Civita tensor. In our approach we consider the J/ip and 
other vector mesons as bound states of constituent quarks and not as gauge fields. 



V. NUMERICAL RESULTS AND DISCUSSIONS 



Some comments about and comparisons of the t dependence of the form factors are in order. The behavior 
of Fpyp(t) = F^lj^it) (in the literature, for t = m 2 D ,, it is called <?D*D7r) and F^ p (t) = F^ D , D (t) in 
the kinematical region is shown in Figs. 4 and 5. In order to be able to compare with other calculations we 
quote the value of F^^^mjj,) (= ,gD*D7r) which is equal to 22. This value is about 2a larger than the recent 
experimental results from CLEO, go'Dir = 17.9 ± 0.3 ± 1.9 [33]. A very small value for 3d*Dtt was predicted 
by the light cone QCD sum rules approach [34]. 

For the Fy^ D ,^(t) form factor, we cannot go on the mass shell due to the presence of an anomalous threshold. 

The dependence of i 7, j/^Q ffD (s,t) on t for different values of y/s is shown in Fig. 6. One can see that the t 
behavior is rather flat. Moreover, the dependence of -Fj/^,D 7r D( s ^) on V$ at t = is shown in Fig. 7. 

In Fig. 8 we separately plot the contributions coming from the box and resonance diagrams for the process 
J/tp + 7r + — > D + + D a . We display the dependence of the cross sections on the variable y/s for each channel 
in Fig. 9. The total cross section is a sum over all channels, 



CTtot(s) = cr D+w (s) + a D , +w (s) + a D+wl7 (s) + a D , +W o(s). 



(27) 
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FIG. 4: Dependence of ^d^dM (solid line) and -Fj/^ D * D (£) (dashed line) on the invariant variable 't" which is the 
D*-momentum squared. 
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FIG. 5: Comparison of our (solid line) form factors ^dj/VidW anc * 'jmd'dW with those (dashed line) obtained in 
[35]. 

which are plotted in Fig. 9 Note that v D+ -pra = g d'+~d°- ^ e P^ ot fJ tot(s) as a function of \fs in Fig. 10. One 
can see that the maximum is about 2.3 mb at y/s » 4.1 GeV. This is close to the result obtained in [4-6]. 
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FIG. 6: Dependence of the form factor -Fj/^ f, ^ D (s, t) on t for different values of a/s. 
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FIG. 9: The cross sections of the processes J/tp + ty + -> D + + D° , J/ip + ty + -> £>* + + L>° and J/i/> + 7T+ ^ L>* + + D* . 
All three curves are shown together in the last plot. 



APPENDIX: EXPRESSIONS FOR FORM FACTORS 



Here we report the analytical expressions for the form factors involved in the calculation of the amplitudes 
J/ip + 7T+ -> D* + + U° and J/i/j + tt+ -> D* + + IF", respectively. 

J/V> + 7T+ £>* + + D° 



2.5 



2 




3.8 4 4.2 4.4 4.6 4.8 5 

Vs (GeV) 



FIG. 10: The total cross section (dotted line) together with the contributions already plotted in Fig. 9 



F vppv( s i l ) = 9J/i> 9d 9d* 9* J da n ■ exp(-w • z u (s, t)) • Nft PV 

iVyppy = — m q m c w 2 + rriq w 2 (2 d\ d 3 — 2 d 2 d 3 — d 3 ) + ml w 2 (-1 + di — d 2 ) 
+m 2 j ^ w 2 (— 4 di d 2 d 3 — 2 d\ d 2 d\ — 2 di d 2 — 2 d\ d\ d 3 
—d\ d\ — 2d\d 3 — d\ d 2 — d\ 
+4 d\ d 2 d 3 + 2d\d 2 + A d\ d 3 
+2d\ d 2 3 + 2d\-2d\ d 3 -d\) 

+771^, w (—6 o?i d 2 g?3 + 2 g?i d 2 d 3 — 2 di d 2 — 4 di d 2 d 3 — 2 di d 2 — 2 di d 3 
+2 c?i + 2 d 2 d 3 + d\ d 2 + 2 d\ d 3 

+2 d 2 d 3 - 3 d 2 d 3 + d 2 + 4 d 2 d 3 - 2 d 2 d 3 + 2 d 2 + 2 d\ d 3 + d\ - d 3 ) 

+m 2 D w 2 (-2 d l d\ + d\ d 3 + 2 d\ d\ - 2 d 2 d 3 

-2 d 2 d 2 + 2 d 2 d\-d 2 d 3 -2 d\ d\ - d 3 - d\ + dj) 

+u w 2 (3 d 1 d 2 d 3 + 2d 1 d 2 + 2 di d\ d 3 + d 1 dj + d 1 d 3 -2 d\ d 2 d 3 -d\d 2 -2 d\ d 3 ) 

+s w 2 (di d 2 d 3 + 2 di d 2 d 3 + d x d 3 + d x d 2 3 - d\d 3 - 2 d\ d\) 

+t w 2 (—di d 2 d 3 — 2 di d 2 d 3 — di d 3 — 2 di d 3 + 2 d 2 d 3 

+3 d 2 d 3 + d 2 d 3 + 2 d 2 d 3 + d 3 + d 3 ) 

+tu (—3 — 8 v di d 3 — 3 v di + 8 v d 2 d 3 + 3 v d 2 + 5 t> d 3 

+3 w + 8 di d 3 + 3 di - 8 d 2 d 3 - 3 d 2 - 5 d 3 ), 

^vppv = +mgW 2 (-2did 2 + d 2 + 2d2) 

+m 2 ji^ w 2 (d l d 2 d 3 + 2 di d| d 3 + 3 di d\ + 2 d l d\ - 2 d\ 

d 2 d 3 -Zd\d 2 - \d\d\ + 2 dl d 2 ) 

+m 2 D w 2 (-2 di d 2 d 3 + 2 di d 2 - 2 di d 2 d 3 + 7 di d\ 
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+4 di d\ - 2 d\ d 2 - 2 d\ d\ + d 2 d 3 - d 2 + 3 d\ d 3 - 4 d\ + 2 dl d 3 - 5 d\ - 2 

+m^. w 2 (—di d 2 d 3 + 2 d\ d 2 d 3 — 2d\d 2 — 2 d\ d 2 d 3 

+3 d 2 d 3 - d 2 d\ + 3 d 2 2 d 3 - 2 d\ d\ + 2 ^ d 3 + d 3 ) 

+u w 2 (-di d 2 - ?> di d% - 2 did\ + 2 d\d 2 + 2 d\d\) 

+s w 2 (—di d 2 d 3 + di d 2 — 2 di d\ d 3 + 2 d\ d 2 d 3 ) 

+t w 2 (2 di d 2 d 3 + did 2 + 2 d 1 d\ d 3 

-di d 3 -3d 2 d 3 -2dl d 3 ) 

+w(l + 8vdid 2 + 2vd 1 -7vd 2 -8vdl-v-8did 2 -2di +7d 2 + 8d?,), 

N VPPV = +m q m c w 2 (2d 3 ) + m 2 q w 2 (-2d 2 d 3 + 2d 2 3 ) + m 2 c w 2 (-l-d 2 ) 
+m 2 ^ w 2 (—di d 2 d 3 - di d 2 - 2 di d\ d 3 - di d\ + di d 3 

+4 di dl + 2 di dl - di + 2 d\ d 2 d 3 + d\ d 2 - 2 d\ d 3 - 2 d\ d\ + d\) 

+m 2 D w 2 (2 di d 2 d 3 — did 2 — 2 di d 2 d 3 

—di d\ + di d 3 + 2 di d 2 — d 2 d 3 — 4 d 2 d 2 

+2 d 2 dl + dl d 3 - 4 d\ d\ + d\ + 2 d\ d 3 + d? 2 + 2 dl) 

+ml w 2 (-di d 2 d 3 -2 di d 2 d^ + 2 di d^ + 2 di dl - d 2 d 3 

+d 2 dj + 4 d 2 dl - dl d 3 - 2 d 2 2 d 2 3 + d 2 3 - 2 d\) 

+u w 2 (—2 di d 2 dl + di d 2 + 2 di dl d 3 + di dl — di d 3 — 2 di dl) 

+s w 2 (di d 2 d 3 + 2 di d 2 dl - 2 d x dl - 2 di df) 

+t w 2 (d 2 d 3 - 2 d 2 d\ + d\ d 3 + 2 d\ d 2 3 - 2 d\ - 2 d 3 3 ) 

+w (-1 + 8 v d 2 d 3 + 3 v d 2 - &v d 3 - 8 v dl + v - 8 d 2 d 3 - 3 d 2 + 4 d 3 + 8 d%), 

-Nyppv = +m q m c w 2 (-2d 2 ) + m 2 q w 2 (-2d 2 d 3 + 2dl)+m 2 c w 2 (l + d 2 ) 
+m 2 /, lll w 2 (-3 di d 2 d 3 -2 di d 2 dl - di d 2 + di dl 

+2 di dl + di d 3 + di + 2 d\ d 2 d 3 + d\d 2 -2 d\ d 2 2 - d\) 

+m 2 D w 2 (-2 di d 2 d 3 -2 di dl d 3 + d x dl + 2 di dl + d 2 d 3 - 2 d 2 dl + 5 dl d 3 

-2 d 2 2 dl - dl + 4 dl d 3 - 3 dl - 2 di) 

+m 2 w 2 (—di d 2 d 3 — 2 did 2 dl — di d 2 + 2 di dl d 3 + di d 3 
+d 2 d 3 - d 2 dl + 2 d 2 dl + d\ d 3 - 4 d\ d 2 z + 2 d\ d 3 - d\) 
+u w 2 (2 di d 2 d 3 + 2 di dl d 3 -d x d 2 2 -2 di dl) 
+s w 2 (di d 2 d 3 + 2 di d 2 dl + di d 2 — 2 di dl d 3 — di d 3 ) 
+t w 2 (d 2 d 3 + 2d 2 d 2 3 - d 2 2 d 3 + 2d 2 2 d 2 3 -2 d\ d 3 ) 

+w(l + 8vd 2 d 3 -vd 2 -8vdl + 2vd 3 -v-8d 2 d 3 + d 2 + 8dl-2d 3 ), 

^VPPV = +m q m c m 2 J/ ^w 2 (-1/2) + m q m c mlw 2 (-1/2) 
+m q m c s w 2 (1/2) + m 2 m 2 w 2 (1) 

+m 2 q m 2 J/4 , w 2 (di d 2 + did 3 + di-d\- 1/2 d 2 - 1/2 d 3 ) 
+m 2 m D w 2 (di d 2 + d 2 d 3 — 1 /2 d 2 — d 2 ) 
+m 2 q ml w 2 (di d 3 + d 2 d 3 - 1/2 d 3 - dl) 
+m 2 q u w 2 (-di d 2 + 1/2 d 2 ) 
+m 2 q sw 2 (-di d 3 + 1/2 d 3 ) 
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+m 2 tw 2 (-d 2 d 3 ) 
+m 2 w (1 — v) 

+m 2 c m 2 J/4 , w 2 (di d 2 + di d 3 + 3/2 di - d\) 

+m 2 c m 2 D w 2 (-1/2 + di d 2 + di + d 2 d 3 - 3/2 d 2 - d^ + d 3 ) 

w 2 (1/2 + di d 3 + 1/2 di + d 2 d 3 + 1/2 d 2 - d\) 
+m 2 uw 2 ( — d\ d 2 — di) 
+m^sw 2 (-di d 3 - 1/2 di) 
+m 2 1 w 2 (-1/2 - d 2 d 3 - 1/2 d 2 - d 3 ) 
+m;;w(2 - 2w) 

+m j m 2 D w 2 (di d 2 d 3 + 2 di d 2 d 3 — 4 di d 2 — 5 d\ d\ 
-2 di d\ + di d^ - 1/2 di + 11/2 d? d 2 + 4d? d^ + 3/2 d? 
-2 dj d 2 - d? - 1/2 d 2 d^ + d 2 + 3/2 d\ + 1/2 d^ - 1/2 d\) 
+m 2 j /ip ml w 2 (2 di d 2 d 3 + 2 di d| d 3 + 1/2 di d 2 - 7/2 di d 3 

-2 di dl + 1/2 di + 7/2 d? d 3 + 4d? d 2 3 - 2 df d 3 - 1/2 d\ - d 2 d 3 - 1/2 d\ d 3 - d 3 + 1/2 d^) 
+m^ /j/) u w 2 (di d 2 d 3 + 3/2 d 1 d 2 + di d| + 1/2 di d 3 -2d\d 2 d 3 

-Ad\d 2 -2d 2 l d 2 2 -d\d 3 -d 2 l +2 dl d 2 + d\) 
+m 2 j^ s w 2 (di d 2 d 3 + 1/2 di d 2 + 3/2 di d 3 + d± d 3 

-2 df d 2 d 3 - 1/2 d? d 2 - 7/2 df d 3 - 2 d\ d\ - 1/2 df + 2 dj d 3 + 1/2 d?) 
+m 2 j /4> t w 2 (-4 di d 2 d 3 - 2 di d 2 d 3 - 1 /2 di d 2 - 2 di d 2 d 3 

-1/2 di d^ -3/2di d 3 -di d 2 3 - l/2d 1 + 2d 2 1 d 2 d 3 

+1/2 d\ d 2 + d\ d 3 + 1/2 d? + d 2 d 3 + 1/2 d 2 d 2 + 1/2 d\ d 3 + 1/2 d 3 + 1/2 d\) 
+m 2 J/ ^w (-3/2 - hvd x d 2 - 5vdi d 3 - 13/2 t>di + 5wdi 

+3/2 wd 2 + 3/2wd 3 + 3/2w + 5di d 2 + 5did 3 + 13/2 di - 5d 2 - 3/2 d 2 - 3/2 d 3 ) 

+ to j/i/j w2 ( — ^1 ^2 d 3 — di d 2 — 1/2 di d^ — d± d 3 — 1/2 di d 3 

+2 d? d 2 d 3 + 3 d^ d 2 + d\d 2 2 + , i d\ d 3 + d^ d^ 

+3/2df-2dl d 2 -2d\ d 3 -b/2d\ + d\) 

+m 2 D m 2 w 2 (—di d 2 d 3 — d\ d 3 + 2 d 2 d 2 d 3 + 1/2 dj d 2 

+d? d 3 + 5/2 d 2 d^ - 2 d 2 d^ + 1 /2 d 2 - d\ d 3 + 4 d\ d\ 

-2d\ d 3 - 1/2 dl - 1/2 d 3 - djj) 

+m D u w 2 (—2 di d 2 d 3 + 3/2 di d 2 — 2 di d^ d 3 

+7/2 di d| + 2 di d^ - 2 d? d 2 - 2 d\ dj + 1/2 d 2 d 3 - 1/2 d 2 

+1/2 d^d 3 -d^ - 1/2 d^) 

+raf, s w 2 (2 di d 2 d 3 - 2 di d 2 d 3 + 1 /2 di d 2 + 2 di d 2 d 3 
+1/2 di d\ + d 1 d 3 -d 1 d\-2 d 2 d 2 d 3 - 1/2 d 2 d 2 - d\ d 3 
-1/2 d 2 d 3 + 1/2 d 2 d\ - 1/2 d 2 - l/2d 2 d 3 - l/2d 2 + 1/2 d 3 ) 
+m 2 D t w 2 (-2 di d 2 d 3 - 1 /2 d Y d 2 -2 d\ d\ d 3 
-1/2 di di + 1/2 d 2 d 3 - 2 d 2 d^ + 5/2 d\d 3 -2 d\ d\ 
+l/2d 2 2 + 2dl d 3 + 1/2 d|) 

+m 2 D w (— 1 — 5vdid 2 — 2vd\ — 5vd 2 d 3 + 4vd 2 
+5vd 2 \-2vd 3 + v + 5d 1 d 2 + 2d 1 + 5d 2 d 3 -4:d 2 -5dl + 2d 3 ) 
+m% w 2 (2 di d 2 d 3 - di d 2 + 2 di d 2 d 3 - 3 di d\ - 2 di d^ + djd 2 + d\ d\ 
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-d 2 d 3 + d 2 d\ + 1/2 d 2 - 3 d^ d 3 + d 2 2 dl + 3/2 d^ - 2 d^ d 3 + 2 d^ + d|) 

+m 2 u w 2 (2 rfi d 2 d 2 - 2 di d 2 d 3 - 1 /2 d\ d\ 

+1/2 di d 3 + di d 2 - 2 <if d 2 d 3 - 1/2 d 2 d 2 - d\ d 3 + l/2d 2 d 3 

-1/2 d 2 d 3 + l/2d 2 d 3 + 1/2 d 3 ) 

+m 2 s w 2 (— di d 2 d 3 — 2 di d 2 d 2 + 1/2 d\ d 3 

+3/2 di d 2 . + 2 did^ - d\ d 3 -2d\d% + l/2d 2 d 3 + 1/2 d 2 d\ + 1/2 d 3 - 1/2 4) 

+m 2 t w 2 (-di d 2 d 3 - 2 di d 2 d 3 - 1 /2 di d 3 - di d 2 

-1/2 d 2 d 3 + 2 d 2 dl~d 2 2 d 3 -2 d 2 d\ + 1/2 d 3 + d 2 + djj) 

+m 2 w (-1/2 - 5 vd 1 d 3 - 3/2 wdi - 5 vd 2 d 3 - 3/2 vd 2 

+5/2vd 3 + 5vdj + l/2v + 5d 1 d 3 

+3/2 di + 5d 2 d 3 + 3/2 d 2 - 5/2 d 3 - 5d§) 

+m* w 2 (di d 2 d 3 + 2 di d 2 d 3 - di d^ - 2 di d| + 1/2 d\ d 3 

+d\ d\ - d 2 d 3 - 2 d 2 dl 1 + 1/2 d 2 . d 3 + d 2 . d 2 . - 1/2 d 3 - 1/2 d 2 . + 1/2 d\ + d 3 ) 

+usw 2 (-di d 2 d 3 - l/2did 2 - 1/2 di d 3 + 2d\ d 2 d 3 + 1/2 d\d 2 + d 2 d 3 ) 

+?itw 2 (2di d 2 d 3 + l/2did 2 + 2di d\ d 3 + l/2d 1 d%- 1/2 d 2 d 3 - 1/2 d 2 d 3 ) 

+uw (1/2 + 5vd 1 d 2 + 2vd 1 -3/2vd 2 -l/2v-5d 1 d 2 -2d 1 + 3/2 d 2 ) 

+u 2 w 2 (-1/2 di d 2 - 1/2 di d 2 2 + d\ d 2 + d\d%) 

+stw 2 (di d 2 d 3 + 2di d 2 d\ + 1/2 di d 3 + di d 2 . - 1/2 d 2 d 3 - 1/2 d 2 d^ - 1/2 d 3 - 1/2 d\) 
+sw{l + bvd 1 d 3 + i/2vd 1 - 3/2 v d 3 - v - 5 di d 3 - 3/2 di + 3/2 d 3 ) 
+s 2 w 2 (- 1/2 di d 3 - 1 /2 di d 2 3 + 1 /2 d 2 d 3 + d 2 d 3 ) 

+t w (-1/2 + 5 v d 2 d 3 + 3/2 v d 2 + 2 1; d 3 + 1/2 w - 5 d 2 d 3 - 3/2 d 2 - 2 d 3 ) 
+i 2 w 2 (1/2 d 2 d 3 + d 2 dj + 1 /2 d 2 . d 3 + d 2 2 d 2 5 ) 
+3-6w + 3w 2 ; 



F PVp(t) = 9TrgD*gD J doA ■ cxp(-wz A (t)) • [wm q m c (-1) + w m 2 (-61) 

+wm 2 (6 2 6 2 - 6 2 + 6?) + ait (61 b\ + 6 2 6 2 - b 2 + &|) 
+wm 2 3 (-bj b 2 ) + v(l + 3 61) - 1 -3 61], 

F pvp(t) = 9.J/->pgD* gD Jdc7 A ■ cxp(-wz A (i)) • [wm 2 (1 - b 2 ) 

+w ml (-3 61 6 2 + h bj + 61 + bj b 2 -b\ + b 2 - b\) 
+w t (-61 b 2 + b 1 bl + b 2 -2bj + b 3 2 ) 
+w m 2 D (61 b 2 - bi b\ - b 2 + 6 2 ) 
+u(-2 + 3 6 2 ) + 2-36 2 ], 

f pvv( u ) = 9t g 2 D* Jda A ■ cxp(-wz A (u)) ■ [w (m q (b 2 - 1) - m c b 2 )], 
f vpv( u ) = 9j/4,gD9D' daA-exp(-wzA{u))-[w(m c (b 2 -l)-m q b 2 )], 



f pvp( u ) = 9* gD* gD J da a ■ cxp(-wz A (u)) ■ [w ((-m c m q - m 2 h) 
+m 2 b\ (62 - 1 + &i) 



+m 2 D b 2 Oi b 2 + bj - 1 + b 2 ) 
-ub 2 1 b 2 )-(l-v){l + 3b 1 )}, 

F vpp( u ) = 9j/^9d Jd°A ■ cxp(-wz A (u)) • [w(-2m c m q b 2 + m 2 (-l + b 2 ) 

+m 2 ^ h {bl-l + h b 2 + h) 
+m 2 D b 2 (b 1 -l + b 2 )-ub 1 b%) 
-(!-«)(! + 36a)] 



J/lj) + 7T+ -> £)* + + D*° 



-^vvpv( s ' ^) 



5i> 5x y cton • cxp(-w • z n (s, i)) • iV^ PV 



(i) 



VVPV 



(2) 



VVPV 



(3) 



VVPV — 
JV VVPV — 
JV VVPV — 



(6) 



VVPV 

iv VVPV 

(8) 

VVPV 
(9) 



N, 



N, 



VVPV 
(10) 



VVPV 



(ii) _ 



VVPV 
(12) _ 



-^VVPV 
iv VVPV 

7V (14) 

JV VVPV 

iV (15) 

JV VVPV 



(16) 



VVPV 



— w m q d 2 , 
w 2 m c d\, 

w 2 m c (1 - d\ + d 2 ), 
w 2 m c (1 — 2d\ + d 2 ), 
—w 2 m c di, 

w 2 (m c (1 - di + d 2 ) - m q d 3 ) , 
-w 2 m q d 3 , 

w 2 dz (m q (2 di — 1) — 2 m c d\) , 

u> 2 <i 2 (m q (1 — 2 c?i) + 2 m c di) , 

w 2 (m q d 2 — m c d\ ) , 

w 2 (ra g d 2 + m c (di — d 2 — 1)) , 

—w 2 m c di, 

w 2 m c (di — d 2 — 1), 

2 u> 2 d 3 (m 9 d 2 - m c (1 + d 2 )) , 

- 2 w 2 d 2 ('Tig c?2 — m c (1 + <i 2 )) 



w 2 m c 



-m 2 q + m 2 ^ 



di [ c?i — - — d 2 — d 3 



1 3 

+ m £>* ( 2 ~~ dlrf 2 ~ rf l ~ ^2^3 + ^ d 2 + d 2 - £^3 



+ m i I - 7: - tW3 - 77 d l - d 2 d 3 - 77 d 2 + d 3 



+ udi{l + d 2 ) + {pi + sdi Q + +* Q(l + d2) + d 3 (l + d2) 

— 2m c w(l — v), 



(17) 



VVPV 



j- m 2 + m 2 7/v> ^-dicfe - cM 3 - di + d 2 + 7j ( d 2 + efe)^ 











1 d 2 - d\ — d 3 H 


-k) 


+ m 2 d 3 I d 3 - di - d 2 H 


-I) 



m 2 D » d 2 



+ u d 2 I d\ - ^ I + s d 3 I d\ - ^ I + 1 d 2 d 3 



2 / 

- w(l-v) (m q + m c ). 



Fpvvit) = 9^9d- J dcr A ■ exp(-wz A (t)) ■ [w(m q (b 2 - 1) - m c b 2 )], 

F-vvvtf) = 9j/^9d* JdcrA ■ exp(-wz A (t)) ■ [4wh b 2 (b x + b 2 - 1)], 

= ffj/v- 3d« / da A ■ exp(-wz A (t)) ■ [2 w b 2 (h - 2 h b 2 - 1 + 3 b 2 - 2 &£)], 

^vvv(*) = ffj/v ffl>* /^a • cxp(-wz A (tj) ■ [w (-m c (m 9 + m c &i) + m} /v , b\ (b x + b 2 - 1) 

+ t 6 2 (6? + &1&2 + &2 - 1) - m% b\ b 2 ) - (1 - u)(6i + 1)], 

^vvv(*) = 9.7/4, 9d- Jdcr A ■ cxp(-wz A (tj) ■ [w {m 2 c (d2 - 1) + m 2 J/xjj (3 b 1 b 2 -b 1 bl~b 1 - b\ b 2 + b\- b 2 

+ t (b, b 2 -b 1 b 2 -b 2 + 2b 2 - b\) - m 2 D (6i b 2 + b 1 bj + b 2 - b\)) + (l-v) (b 2 - 1)], 

^vvv(*) = S'J/V 3r>* • exp(-u7z A (i)) ■ [w (m c (m g - m c &i - m c b 2 + m c ) 

+ m2 //v fe i (~ 2 fe 2 + b 2 2 + 1 + 2 6i 6 2 - 2 6i + 6?) 

+ i 6 2 (-2 &i + 2 &i b 2 + b\ + 1 - 2 fo 2 + 6^) + m 2 D h b 2 (2 - b 2 - h)) + (1 - v) (2 - h - b 2 )], 

Fyvv(t) = 9j/4,9d* Jda A ■ cxp(-wz A (t)) ■ [w (m c (2m q b 2 - m c b 2 + m c ) + m 2 ^ b\ {-b\ + 1 - b x b 2 - 

+ tb 2 (-61 b 2 + b 2 - b 2 ) + m 2 D 61 b 2 ) + (1 - v) (62 + 1)], 

Fppy(t) = 9-k9d9d* fdcr A ■ exp(-wz A (t)) ■ [w (m q (m c - m q bi - m q b 2 + m q ) 

+ m 2 7/j/) h (1 - 26i - 2fo 2 + b\ + 2hb 2 + b\) 

+ 1 6 2 (1 - 26i - 26 2 + bj + 2b x b 2 + bl) + m 2 D h b 2 (2 - h - b 2 )) + (1 - v) (4 - 3 61 - 3 62)] 
^vpv(*) = 9^ / rf °A • cxp(-w7z A (i)) • [to (m c - 6 2 (m c - m 9 ))]. 
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